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Vipin Kerala Varma and Hartmut Monien 

Bethe Center for Theoretical Physics, 
Universitat Bonn, Germany 

(Dated: November 27, 2012) 

Lattice Green's functions (LGF) and density of states (DOS) for non-interacting models on 3 
related lattices are presented. The DOS and LGF at the origin for the kagome and diced lattices 
are rederived. Furthermore, from the form obtained for the DOS of the hyperkagome lattice, it 
seems highly suggestive that a closed form expression can be obtained if a particular complicated 
yet highly symmetrical integral can be solved. The results have been checked numerically. The 
continuum of energy states in the hyperkagome lattice and the appearance of a zero in the kagome 
lattice may be explicitly seen from the formulae obtained. 



Lattice Green's functions (LGF) of systems are 
ubiquitous^ in solid state physics appearing in problems 
of lattice vibrations, of spin wave theory of magnetic sys- 
tems, of localized oscillation modes at lattice defects, and 
in combinatorial problems in lattices^. The LGF G(t, r) 
in the complex energy variable t = s — ie between the 
(arbitrary) lattice origin and the point r is defined as the 
solution to the inhomogcnous difference equation 



[2t-H(r)]G(t,r) = 2S, 8 



(1) 



corresponding to the linear, Hermitian, time-independent 
operator H(r). This operator is generally taken to be the 
discrete Laplacian for lattice problems. Its action on an 
n-variable function 0(r = {r\, ri • • ■ r n }) is given by 



(r + A) -z<}>(f), 



(2) 



where N denotes the set of nearest neighbours of r and 
z = \Af\ is the number of such elements. In subsequent 
analysis, we will neglect the second term in Eq. (2) be- 
cause it merely adds a constant shift to the energy vari- 
able t. With that Eq. (1) defining the LGF is expressed 



[2tG(t,f)- J2 



= 2^,o- 



(3) 



The above equation may be readily solved by taking its 
Fourier transform into k-space, noting that G(t, r+ A) O 
e ( k A )(5(i j k), where G denotes the Fourier transformed 
variable. For a simple cubic lattice in e?-dimensions with 
unit lattice distances, this can be written, after inverse 
Fourier transforming, as 



G(t,r) 



( — Y 



-7T t — ^(k) 



(4) 



where uj(h) = J^. cos (ki). The general form for the LGF 
for other lattices like body-centred cubic, face-centred cu- 
bic remain the samei. In many cases such integrals may 
be exactly solved in terms of the elliptic integrals of the 



first and second kind. Once a solution for the LGF is ob- 
tained, we may obtain information about the eigenvalues 
and eigenfunctions of the operator H(r). And any solu- 
tion of the inhomogenous equation may be constructed 
using G(t, f). An important quantity of the lattice prob- 
lem is the density of states (DOS) p(s), which will be a 
focus of this paper, given as 



p(s)= lim -3[G(t,0)], 

e-s-+0 7T 



(5) 



where S refers to the imaginary part. 
In this brief report, we present results for three related 
lattices at r = 0. In the first section, we derive the 
LGF and DOS for the kagome lattice, which is shown 
in Fig. 1(a). In the second section we perform similar 
calculations for the diced lattice, which is shown in Fig. 
1 (b) . Some of these results were obtained earlier— but we 
repeat it here to correct a minor error in Ref. [4]. And 
in the final section, we attempt to obtain a closed form 
solution for LGF and DOS for the hyperkagome lattice, 
which is shown in Fig. 2(a). We could not solve this 
exactly in a closed form but numerical evaluation of this 
integral to obtain the DOS does, however, compare well 
with a previously published purely numerical evaluation 
of the DOS. 



I. KAGOME LATTICE 

The kagome lattice has 3 atoms per unit cell labelled 
A,B,C in Fig. 1(a). The side-length of the simple tri- 
angle a — 1 and the vertical height b = We use 
the general formalism presented before to calculate the 
LGF at the origin using Eq. (3) and the DOS from the 
imaginary part using Eq. (5). The difference here is that 
we need to consider LGFs at and between the different 
sublattices a,/3 denoted as G a/3 . Then the right hand 
side of Eq. (3) is multiplied by an additional Kronecker 
delta factor S a t p and Eq. (5) is modified to 



(0) 



2 



where N denotes the number of sublatticcs. 
Taking the Fourier transform of Eq. (3) with respect to 
r, the lattice Green's functions, evaluated at the origin, 
for the ^-sublattice with respect to the other 3 sublattices 
may be written as 

t Q A u _ qBv cqs ^ + bk ^ _ qCv cqs ^ _ bk ^ = 5a ^ 

-G Av cos (ak x + bk y ) + tG Bv - G Cv cos (2ak x ) = 5 Bu 
-G Av cos (ak x - bk y ) - G Bv cos (2ak x ) + tG Cu = 8 Cv , 



Pkag.(s) 



3tt L 



1 



where 



2s 2 



t+1 



^[2 



2t- 



(9) 



2s — 1, L)' k = cos (Aak x ) 



2 cos (2ak x ) cos (2bk y ) and T denotes the Fourier trans- 
form of the k variable within the limits ±7r/a, ±7r/&. 
That is, 



Pkag.( s ) 



S(s + 1) 



+ -|2s-l|p tri (r / ). (10) 





where v = A, B, C as before. The LGFs G a ° depend on 
(t, k x , ky) but has been suppressed in the above equation 
for notational simplicity. It is a straightforward task to 
solve for G a/3 from the above set of linear equations by 
setting particular values for v. Then we may write the 
"diagonal" lattice Green functions i.e. a = (3 as 

G AA = [t 2 - cos 2 (2ak x )]/V, 
G BB = [t 2 - cos 2 (ak x - bk y )]/V, 
G cc = [i 2 - cos 2 (ak x + bk y )]/V, 

(8) 

where T> = (t + l)[t(t — 1) — cos (2ak x )(cos (2ak x ) + 
cos(2bk y ))]. The "off-diagonal" lattice Green functions 
may be similarly calculated as well from Eq. (7). 

From the knowledge of the Fourier-space diagonal 
LGFs G aa in Eq. (8), the DOS p kag .(s), as a function of 
the real energy variable s = $t[t], for the kagome lattice 
may be calculated using Eq. (6) 



It may be checked that the spectral density integrates to 
one i.e. J , /9kag.( s ) ds = 1. The expression in Eq. (10) 
checks with that in Ref. [5] which provided no derivation 
of the result. We additionally have obtained the diagonal 
LGFs, in Fourier space, for the kagome lattice in Eq. (8). 



II. DICED LATTICE 

The diced lattice is the dual to the kagome lattice 
and is shown in Fig. 1(b) with its three sublattices again 
labelled as A, B, C . The length of segment AB is taken 
to be one. Here we merely recapitulate the work of Ref. 
[4] in slightly different form, to correct a very minor error 
in the density of states. The calculations proceed as in 



(a) 



(b) 



FIG. 1: (a) Diced lattice, (b) Dual to the diced, kagome 
lattice. 



the previous section and the missing details may readily 
be worked out or obtained from Ref. [4]. The LGF at 
the origin G a,a (t,k x ,k y ) - where a = (A,B,C), (k x ,k y ) 
are the Fourier space components as in Section I - for the 
3 sublattices A,B,C are given by 



G (t, k Xl ky) 
G (t, k x , ky) 
G {t, k x ,ky) 



4f 



4t 2 -(l+ 7 2 K' 

4* 2 - 7 2 ^fc 
t[At 2 - (1 + 7 2 K] ' 

4t 2 - LO k 
£[4t2-(l + 7 > fc ]; 



(11) 



where u>k = 2cos(2ak x ) + 4 cos (ak x ) cos (bk y ), 7 is the 
A <-> C hopping amplitude in units of the A <-> B hopping 
amplitude, which is taken to be 1. a = V3/2 and 6=1/2 
are lattice distances in units of the inter-atomic distance. 
Transforming to real space i.e. G — >• G, the DOS as a 
function of the real energy variable s = 3l[t] is given by 



Pdiced(s) 



3^E G ° 



8(8) 



+ 



4s 



3tt(1+7 2 ) 



9f[Gtri(r > 0,0)] > (12) 



where r 



4s J -3(l+7 2 ) 



2 (i +1 2) — an d Gtri(t, x, y) is the LGF for the 
triangular lattice; the latter is known from other worls^. 
Then Eq. (12) simplifies to 



PdicecH 5 ] — 



4s 



3(1 + 7 2 ) 



Ptri(r). 



(13) 



The DOS for the triangular lattice derived and expressed 
in Ref. [5] in terms of the complete elliptic functions of 
the first kind K(m), with m being the parameter, may 
be rewritten, in the whole range —1.5 < r < 3, as 



Ptri(r) 



5R[A^l-l/fc(r) 2 )] 



where k{x) 



tt 2 (2t + 3)3 



3 



Here too, it may be checked, say for a particular 

3 

value of 7 = 1, that Pdiced(s) ds = 1. In Rcf. [4], 

the Dirac-dclta term was not taken into account and the 
spectral density did not integrate to 1; it is clear that 
the flat band in the diced lattice must give rise to such a 
term. 



uated; a rationale for this proposition was later giveni. 
Glasser's canonical form for the Green's function is 



where 



G — — r 

7T J 



dxdydz 



o t — F(cosx, cosy, cos 2) 



(14) 



III. HYPERKAGOME LATTICE 

The hyperkagome lattice can be considered to be a 
three dimensional version of the kagome lattice. Here, 
there are corner-sharing tetraheda with one of the atoms 
in each tetrahedron removed. The unit cell, with 12 
atoms, is cubic and is shown in Fig. 2(a). Each side of the 
cube is taken to have unit length. For example, the lattice 
point labelled 3 has (x,y,z) coordinates (3/4,1/4,1/2), 
with respect to the origin O, with the nearest neighbours 
labelled 1, 2, 4, 5. 

There currently is hardly any exactly known result for 
this lattice, even for non-interacting models. However 
it was conjectured^ that lattices with cubic cells should 
have a lattice Green's functions expressable in a "canoni- 
cal" integral form, which, in general, can be exactly eval- 



F(x,y, z) = aix + a2y + a3Z + a4xy + a5yz + aexz + ajxyz. 

Here, based on calculations for the hyperkagomc's 
LGF and DOS, we obtain a form that resembles the 
canonical form of Glasscr. This we solve numerically and 
compare with other purely numerical results. 

Corresponding to each of the 12 lattice points, each 
of which have 4 nearest neighbours, there will be one 
LGF equation. These may be obtained in the same way 
that led Eq. (3) to Eq. (7). The system of equations for 
the a sublatticc, using the corresponding LGF column 
vector 0f 2xl = [G 1,Q • • • G 12 ' a ] T 7 is given by 



CQ a = 2[6i 



■ £l2,a]l2xl! 



(15) 



where the 12x12 coefficient matrix C is given by 



C = 



I 2t 


£, — x — z 
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2t 


£,x+y 

















S— J/— J8 
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2t 
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Zx — z 
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2t 


£x-y 








x — z 
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and the symbol £,±^±v = exp^ ±fcf ' ±fc "^. Now there arc Here the inverse Fourier transform T 1 , at r = 0, of a 



12 systems (for each value of a) of 12 linear equations, 
which is cumbersome to solve by hand as done in previous 
2 sections for smaller systems. To obtain the diagonal 
LGFs, we use Mathematical to obtain each G aa . With 
these solutions we can use Eq. (6) to express the DOS 
Phyporkag. for the hyperkagome as 



function g is given by J-~ 1 g = -^s JJJq <?(k)dk. The a* 
and bi parameters are only functions of t with no depen- 
dence on the lattice momenta i.e. 



_ 1 ^[r-lf ^" ' 81 '" 2 '"^ 1 1 
Phypcrkag. - { p^ ^ ^ I f + ^ 



(16) 
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FIG. 2: (a) Unit cell, from Ref. [11] with authors' permission, of the hyperkagome lattice with blue lines and filled circles; the 
origin is marked with a red O. The additional empty circles indicate the underlying pyrochlore lattice, (b) Density of states 
for the hyperkagome, from evaluating Eq. (16), compared with numerical results from Ref. [9]. Note that the latter had to be 
multiplied by a factor of 2 to account for their doubling of the energy scales compared to ours. 



a = 3(-l + 12t + 72t 2 - lUt 3 - 256t 4 + 448i 5 + 128i 6 - 384t 7 + 128i 8 ), 
ai = -4(1 + 7t - 6t 2 - 20t 3 + 16i 4 ), 
a 2 = 4(1 + 3t - 6t 2 ), 
a 3 = 12£(1 - 2t), 

b = 1 - 16* + H2i 2 + Mt 3 - 544i 4 + 256i 5 + 512t 6 - 512i 7 + 128i 8 , 
h = -16t(l + t- U 2 + 2t 3 ), 
b 2 = — 16*(t — 1), 

6 3 =4(-l+4t-4t 2 ). (17) 
And the P function is given by 

33 33 
P(a, b,c,d) = a + & / J cos {xi) + c \_] cos (xi) cos (xj+i) + dTJ cos (xi) + 2_. cos (2xj), 



Numerical evaluation of the hyperkagome's DOS^ is com- 
pared with a numerical evaluation of the integral in Eq. 
(16) and is shown in Fig. |2(b)| Note that Eq. (16) cor- 
rectly shows the appearance of a flat band at t = — 1 as 
indicated by the simple pole. And as noted from numer- 
ical calculations in Ref. [9], there is a drop to zero in 
the kagome lattice's DOS for the non-interacting Hub- 
bard model (which is merely the second-quantized form 
of the operator H(r) ) but a continuum of energy states 
in the hyperkagome. We may see this explicitly from Eq. 
(9) with the appearance of a zero at 2t — 1 = 0; no such 
apparent zero occurs for the case of the hyperkagome in 
Eq. (16), verified by numerical integration of the same as 
shown in Fig. |2(b)| A continuum of states is, of course, 



indicative of conductive behaviour. 

Many similar integrals - which can be exactly solved - 
appear in the literature^^ for the lattices with a cubic 
unit cell but not with all the sums and products of H 
appearing together. In particular, the last summations 
of cos 2x in the P-function do not appear in the canon- 
ical form in Eq. (15), (16). It is unclear if these terms 
might suggest a generalization of the canonical form or 
if integrability is destroyed by their presence. We have 
currently found no way to exactly solve the Eq. (16). 

To summarize, we have calculated lattice Green's 
functions and density of states for three related lat- 
tices: the kagome, diced and hyperkagome. For the 
hyperkagome lattice, however, we are able to reduce the 



5 



obtained integral to a form approximately like Glasser's 
canonical form but with some additional terms that 
appear to destroy intcgrability. Due to the relatively 
succinct functional form obtained for the integral and 
the numerical correspondence of this result with other 
purely numerical evaluations for the density of states, 
we believe that it might be plausible to obtain a closed 
form solution. Exact expressions for lattice density of 
states ought to be useful in dynamical mean field theory 
calculations^. 
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